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We present a detailed study of nonequilibrium Joseph- 
son currents and conductance in ballistic multiterminal SNS- 
devices. Nonequilibrium is created by means of quasiparticle 
injection from a normal reservoir connected to the normal part 
of the junction. By applying a voltage at the normal reservoir 
the Josephson current can be suppressed or the direction of 
the current can be reversed. For a junction longer than the 
thermal length, L 3> £t, the nonequilibrium current increases 
linearly with applied voltage, saturating at a value equal to 
the equilibrium current of a short junction. The conductance 
exhibits a finite bias anomaly around eV ~ Tivf/L. For sym- 
metric injection, the conductance oscillates 27r-periodically 
with the phase difference <f> between the superconductors, with 
position of the minimum (0 = or n) dependent on applied 
voltage and temperature. For asymmetric injection, both the 
nonequilibrium Josephson current and the conductance be- 
comes 7r-periodic in phase difference. Inclusion of barriers at 
the NS-interfaces gives rise to a resonant behavior of the to- 
tal Josephson current with respect to junction length with 
a period ~ Xp. Both three and four terminal junctions are 
studied. 

PACS: 74.50.+r, 74.20.Fg, 74.80.Fp 



I. INTRODUCTION 

The art of controlling Josephson current transport 
through mesoscopic superconducting junctions poses 
many challenges for theory and experiment from both 
fundamental and applied points of viewEI Control of 
Josephson current requires multi-terminal devices - su- 
perconducting transistors. One example is the Joseph- 
son field effect transistor (JOFET),at3 where control of 
the Josephson equilibrium current is imposed via an elec- 
trostatic gate. Another solution is to connect the nor- 
mal region to a normal voltage biased reservoir. Recent 
progress in fabrication of superconducting junctions has 
brought forward a number of interesting multiterminal, 
structures, e.g. 2DEG,.iunctionsuT3 metallic junctionsJj'El 
and high-Tc j unctions .0 

Injection of electrons and holes allows nonequilibrium 
quasiparticle distributions to be maintained in the N- 
rcgion, making it possible to control the nonequilibrium 
Josephson current. The problem of nonequilibrium cur- 
rent injection in ballistic junctions is of particular inter- 
est: the Josephson, current is transported through bands 
of Andreev levels This provides means for achieving a 
dramatic variation of the Josephson current. 



The purpose of this paper is to provide a broad de- 
scription of Josephson current transport through ballis- 
tic SNS junctions under conditions of nonequilibrium in 
the normal region due to contact with a voltage biased 
normal reservoir. Connection of the normal part of an 
SNS junction to a normal electron reservoir gives rise to 
broadening of the Andreev bound levels. Van Wees et 
al.E£l were the first ones to consider this broadening in 
perfect SNS junctions and to describe essential aspects 
of the variation of the Josephson current with voltage, 
in terms of nonequilibrium popu.l3t.ian of Andreev levels. 
Moreover, Wendin and ShumeikcElrEj have predicted that 
nonequilibrium filling of Andreev levels may reveal very 
large Josesphson currents with different directions, and 
that pumping between levels could reverse the direction 
of the Josephson current. This has been investigated in 
detail by Bagwell and coworkers^ and by Samuelsson et 
alH 

Recently, Morpurgo et al.lia observed Andreev levels by 
using the injection lead as a spectroscopic probe. Sup- 
pression of the Josephson current due to injection has 
been demonstrated in both ballisticO and diffusive SNS 
junctionsu. The physical mechanism of the effect in dif- 
fusive junctions is essentially the same as in ballistic 
junctions£3 Very recently, Baselmans et. al.Ll were able 
also to reverse the direction of the Josephson current. _ 

A decisive step beyond the work of van Wees et al.EHl 
was taken by Samuelsson et alii-l, who showed that an 
essential aspect is the ability of the scatterer at the in- 
jection point to shift the phases of the quasiparticles. 
In such a case, the connection to the injection lead also 
affects the form of the wave function of the Andreev res- 
onances, and therefore affects Josephson currents flow- 
ing through the resonances. As a result, modification of 
the Josephson current under injection does not reduce 
to the effect of non-equilibrium population. This is par- 
ticularly dramatic for long junctions, where the equilib- 
rium Josephson current is exponentially small at finite 
temperatureO 

In contrast, this anomalous nonequilibrium Josephson 
current does not depend on the length of the junction 
(long-range Josephson effect). This means that, in prin- 
ciple, a dissipationless current of the order of the equi- 
librium Josephson current of a short junction can be re- 
stored under conditions of filling up all the Andreev levels 
in the gap. The effect is most pronounced in junctions 
with a small number of transport modes. This opens 
up the possibility for a new kind of Josephson transis- 
tor where the supercurrent is turned on when the gate 
voltage is switched from eV = to eV = A.E3 
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The complete picture of the nonequilibrium current 
also includes the current injected into the junction. This 
injection current is dependent on the properties of the 
Andreev levels in the junction. It therefore provides in- 
formation on the nonequilibrium Josephson current. We 
found that it is closely related to the anomalous current 
and has similar properties. The injection current has 
also in itself been at the focus of great interest in recent 
literature O 

The paper is organized as follows. In section II 
we present a general discussion of the currents in a 3- 
terminal SNS device. In section 77/ we describe our 
model based on the stationary BdG equation. We de- 
rive all currents in the case of a three terminal junction 
without barriers at the NS-interface in section IV. In 
section V we discuss the equilibrium and nonequilibrium 
Josephson currents, both in a short and long junctions. 
The effect of barriers at the NS-interfaces is discussed in 
section VI and the injection current and the conductance 
are analyzed in section VII. In section VIII we discuss 
the four terminal junction and how it differs from the 
three terminal one. Finally, in section IX we present our 
conclusions. 



II. NONEQUILIBRIUM JOSEPHSON CURRENTS 

We will consider two junction configurations: 3- and 
4-terminal (see Fig. Q). The normal part of the junc- 
tion is inserted between two superconducting electrodes. 
The superconducting electrodes are connected with each 
other to form a loop and the magnetic flux threading the 
loop allows us to control the phase difference <fi = 4>r — 4>l 
across the junction. 




FIG. 1. A schematic picture of the three terminal SNS 
junction setup under consideration, with a normal reservoir 
attached to the normal part of the junction. The normal reser- 
voir is connected to the superconducting loop (grounded) via 
a voltage source biased at V . The right figure shows a close-up 
of the junction area with the arrows showing the direction of 
the current flow in the junction 

We consider a junction in the ballistic limit, i.e when 
the length L = L2 + L3 of the normal part of the junction 
is shorter then both the elastic and inelastic scattering 
lengths, L -c Z e , Zj. The 3-terminal configuration is an el- 
ementary structure which gives all necessary information 
for understanding also the properties of the 4-terminal 



junction, to be discussed below. We use a simplified de- 
scription of the connection point, modeling it by a scat- 
tering matrix S thatxpnnects ingoing and outgoing wave 
function amplitudesEHl 

*out = SV in , (1) 

with 

S = K/i r d , (2) 
d r ) 

where r and d are reflection and transmission amplitudes 
for scattering between lead 2 to lead 3 and y/e is the 
scattering amplitude from the injection lead 1 to lead 2 
or 3. In a multichannel treatment, r, d and e become 
matrices describing the scattering between the channels. 
In this paper we however choose to consider a single-mode 
structure. 

In the junction presented in Fig. [l], the current I\ in- 
jected into the junction from the normal reservoir splits 
at the connection point. At the NS-interfaces, the nor- 
mal current is converted into a supercurrent. The super- 
current flows around the loop and is drained at a point 
connected to the normal reservoir via a voltage source bi- 
ased at voltage V . There are two major questions about 
the currents: (i) what is the current I\ in injection elec- 
trode 1 as function of the applied voltage, and (ii) how 
is the current split between the ar ms 2 and 3. The first 
problem has been discussed earlier ,EM3 the picture is the 
following: due to Andreev quantization the problem is 
equivalent to a resonant transmission problem. For weak 
coupling to the normal reservoir, e <C 1, the probability of 
an incoming electron to be reflected is large unless its en- 
ergy coincides with an Andreev level. In such a case, the 
electron is back scattered as a hole which produces a cur- 
rent density peak. The current as a function of applied 
voltage between the normal reservoir and the junction 
(IVC) thus increases stepwise, typical for resonant trans- 
port, with position and height of the steps depending on 
the phase difference between the superconductors. 

The current distribution among the left and right arms 
of the junction is also phase dependent. However, there is 
a less trivial aspect of the problem related to the Joseph- 
son current in the loop. There is no possibility to distin- 
guish the Josephson current which flows along the loop 
(as the result of an applied phase difference) from the 
split injection current except in the limit of weak coupling 
to the external reservoir. In the limit e < 1 the injec- 
tion current (~ e) vanishes while the Josephson current 
remains finite. This allows us to separate the problem of 
the Josephson current under injection from the problem 
of splitting of the injection current. 

The scattering states carrying the current can qualita- 
tively be described as electrons or holes entering the SNS 
junction from the injection lead 1, being split at the con- 
nection point, scattered back and forth in the junction by 
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Andreev reflections at the NS-interfaces and normal re- 
flections at the connection point, and then finally leaving 
the junction, having effectively transported current from 
one superconductor to the other. When the lifetime of 
the Andreev resonances is smaller than the inelastic scat- 
tering time in the junction, the quasiparticle distribution 
in the normal region is determined by the Fermi distri- 
bution function of the normal reservoir, and the current 
in the leads j = 2 or 3 from injected quasiparticles can 
be written 



dE^ri 



-\-lj7l ) 



(3) 



where ij is the current density for injected electrons 

(holes) and n et - h ^ = uf(E ± eV) are the Fermi distri- 
bution functions in the normal reservoir, with np 
1 + exp(E/kT)]^ 1 . This current can conveniently be 
rewritten 



dE 



(n' 



ii 



where 



and i 



(4) 

i e — i h . Quasiparticles are 
also injected from the superconductors for energies above 
the superconducting gap. Since the superconductors are 
grounded (V = 0), the current from the superconductors 
is an equilibrium current. This current plus the current 
I + = 1-2 = I3 injected from the normal reservoir in 
absence of applied voltage, is the total equilibrium cur- 
rent. Applying a bias voltage (V 7^ 0), / + becomes the 
nonequilibrium current due to population of the empty 
Andreev levels, giving rise to current jumps when the in- 
jection energy eV equals the Andreev level energies (see 
Fig. ||). This makesjit possible to probe the energy of 
the Andreev levels-EMO 
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FIG. 2. The current voltage characteristics (IVC) for I 
(upper) and I~ (lower) for a junction with seven Andreev 
levels for < E < A. The currents jump every time the 
voltage eV is equal to the energy of an Andreev level, typical 
for resonant transport 



The J - part of the current is entirely nonequilibrium 
current. It partly consists of the injection current; how- 
ever, there is also a component which does not vanish 
in the limit of weak coupling to the .reservoir: we call 
this the anomalous Josephson currents^ This current re- 
sults from a different form of the Andreev resonance wave 
functions in the open junction compared with the wave 
functions of true Andreev bound states. The origin of 
the anomalous current can qualitatively be described by 
considering the lowest order quasiparticle classical paths 
which contribute to the resonances in transparent junc- 
tions (i? <C 1) with perfect NS interfaces. 




E/A 

FIG. 3. The charge current density for two resonant An- 
dreev levels for injected electrons i e (dotted) and holes i h 
(dashed), their sum i + (solid) and difference i~ (dash-dotted). 
Note that the difference current i~ has the same sign for both 
resonances. Inset: Two lowest order paths for an injected 
electron (solid) or a hole (dashed) at a resonance. The grey 
ellipse denotes the effective scatterer due to the three lead 
connection. The difference of the currents due to these pro- 
cesses is proportional to Im(rd*) sin0, the first order term of 
the anomalous current 



Consider a resonant state where the most of the elec- 
trons move to the left and the holes to the right, only 
a fraction of them travellling in the opposite direction 
due to normal scattering at the connection point. An 
injected electron gives rise to a left going electron in lead 
2 with the amplitude l + e !fa d*e" lfe r with 4> = 4>r~4>l 
(not taking electron and hole dephasing and the energy 
dependent phase picked up when Andreev reflecting into 
account) thus giving a contribution to the current of the 
order 1 + RD + Re(rd*e^) (see inset in Fig. ||). Corre- 
spondingly, an injected hole gives rise to a rightgoing hole 
in lead 3 with amplitude 1 + e~ z ^ L de l ^ R r* and a contri- 
bution to the current of order 1 + RD + Re(rd* e~ 1 ^) (see 
right figure in inset in Fig. ^|). The difference current i~ 



thus contains a part proportional to Re[rd* (e 



i(f) _ 



)] 



2Im(r<i*) sin(0), which is the leading term in the anoma- 
lous current. At a resonant state where the particles 
move in the opposite direction, i.e the electrons to the 
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right and the holes to the left, we find from the same argu- 
ments that the anomalous current is again proportional 
to 2Im(rd* ) sin(</>) , with the same sign. The anomalous 
current thus flows in the same direction for all resonances, 
in contrast to the equilibrium Josephsson current which 
changes sign from one level to the next. The IVC for 
I~ is thus a staircase, as shown in Fig. g, saturating 
at eV > A due to the absence of sharp resonances for 
energies above the superconducting gap. This has a dra- 
matic effect on the long range properties of the Josephson 
current. 

For a long junction (L > £ = hvp/A), the IVC in Fig. 
U becomes dense, since there is a large number ~ L/£o of 
Andreev levels in the junction. The spacing between the 
Andreev levels is ~ hvp/ L, so at temperatures exceeding 
the interlevel distance, the current I + is averaged to zero 
while I~ is reduced to a smooth ramp function. We 
thus get a current I~ that increases linearly with voltage 
up to eV = A and saturates at a level of the order of 
the equilibrium Josephson current of a short junction, 
I ~ eA/h. This current is independent of the length of 
the junction, since there is a large number of levels ~ L 
each carrying a current ~ 1/L. 



III. CALCULATION OF THE CURRENT 

A. General formulation 

We consider a three-terminal junction with asymmetric 
current injection {L2 7^ L3) and perfect transmission at 
the NS interfaces. The junction can be described by the 
stationary 1-D Bogoliubov-de Gennes (BdG) equations 



H A 
A* -H 



which gives E^-^s a departure from Ep. We apply the 
approxinmtiorJH-J with A(x) constant in the superconduc- 
tors and zero in the normal region. 



Ae^ L x < -L 2 
A(x) = { -L 2 <x<L 3 
Ae^ R x > L3 



(6) 



where the phase difference between the superconductors 
is 4> — 4>r — 4>l- We can then make an ansatz with plane 
waves in the different regions of the junction. For positive 
energies E > we put in the normal regions j — 1, 2, 3, 



+C 4 



e lk ° x + c]~ 



-ik c 



and in the superconductors j = L, R 



ue t<t> 3 
v 

v 



iq x 



jh,- 



-iq x 



The coherence factors u and v are defined as 



«(+),«(-) = 



s(i±e/g) e> a 

h(E±Q/A E < A 



(8) 



(9) 



where ^ = y/E 2 - A 2 for E > A and g = WA 2 - E 2 for 
E < A. The wavevectors are q e ' h — Jlmjhf ! y 'Ep ± ^ in 



the superconductors and and k e,h — J 2m/h 2 y/Ep ± E 
in the normal regions. The wavefunctions are matched 
at the NS-interfaces and at the injection point. The 
three-terminal .injection point is modeled by the scatter- 
ing matrixEE! given by Eq. (|J). The scattering am- 
plitudes e (0 < e < 0.5), d and r obey the relations 
Re{rd*) = -e/2 and D + R = 1 - e (D = \d\ 2 ,R = \r\ 2 ) 
due to the unitarity of the scattering matrix. Moreover, 
Im(rd*) = cry 1 RD — e 2 /4, with a = ±1 dependent on 
the phase of the scatterer. For simplicity the coupling 
parameter e is chosen real and positive. The scattering 
amplitudes are assumed to be energy independent, which 
gives the scattering matrix for hole wavefunction ampli- 
tudes S h = S*. 

Assuming A <^ Ep we make the approximation q e = 
q h = kp in the superconductors and k e — k h = kp in 
the normal region except in exponentials where we put 
k e - h — kp ± E/(Hvp). At energies E < A, only elec- 
trons and holes from the normal reservoirs are injected 
in the junction. For E > A quasiparticles from the su- 
perconductors are also injected. The current density in 
the three normal regions, which is what is needed to cal- 
culate all currents in the junction are calculated using 
the quantum mechanical formulaEj 



ij(E) = ZQcfr \cj-r |c+'Y + IcjT)- (10) 



,h\2\ 



We now define energy dependent phases 62,3 = 7 — P2.3 
in each of the leads 2 and 3, consisting of the phase 
7 = arccos(£'/A) picked up by the electrons and holes 
when Andreev reflecting, and the dephasing ^2,3 = 
(k e — k h )L2,z = 2EL 2 ^/ (%vp) of the electrons and holes 
while propagating ballistically through the normal re- 
gion. Furthermore, it is convenient to separate out the 
specific features of asymmetry by introducing sum phases 
29 = 8 2 + 3 , (3 = 02 + 03, an d the difference phases 
X = 82 — O3, defining essential phase parameters charac- 
terizing the junction, 

6 = 7 - j9/2 = arccos(£7A) - EL/(hv F ) (11) 
X = fo-P2 = 2El/(hv F ) (12) 

where L = L2 + L3 and I = L3 — L2 
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The current densities of the scattering states in leads 
2 and 3 from electrons i% 3 and holes i% 3 are then given 
by 



■ e.h 



hz 



{2D sin</>sin26» 



± 



a2 v / RD-e 2 /4 sin </>(cos x — cos 29) 



+e [1 - cos(2 7 - /? 2 ) + cos ^(cos x ~ cos 26)}}} (13) 



= — — — {2D sin <j> sin 20 
h Z 



± 



V V4 X — cos 20) 

-e [1 - cos(2 7 - /3 3 ) + cos 0(cos x - cos 26)}}} (14) 

where 

Z = [(1 - e) cos 26 - i? cos x - cos 0] 2 + e 2 sin 2 26> 

(15) 

From Eqs. (|l^) and ([lj]) it follows that the sum of the 
electron and hole current densities, i + = i e +i h , are equal 
in leads 2 and 3, giving the sum current density 



/ ., = 14 = — ^-^{£>sin(?!)sin20}. (16) 

ft Zj 



'3 



The difference current densities i~ = i e — i h in leads 2 
and 3 are not equal, however. We therefore define the 
anomalous current density i a as that part of the differ- 
ence current density which survives in the limit e — > 0, 



RD - e 2 /4sin</)(cosx - cos26»)| (17) 



4e e 

r ~hz 



The injection current density ij n j = i 3 — i 2 is given by, 



4e e 



— — { sin 2 x + (cos X + cos 0) (cos x — cos 26) } 



(18) 



B. Weak coupling limit 

Throughout the paper we will mainly discuss the sit- 
uation when the normal reservoir is weakly coupled to 
the normal part of the junction, e <C 1. In this limit 
the Andreev resonances are very sharp and the current 
densities are calculated by evaluating the expression e/Z 
appearing in the Eqs. (JT^)- ( 18), in the limit e — > 0. This 
is done in detail in Appendix B, and gives [see Eq. ( |B6] )1 



Hm - = y 

Z ^ 



n.± 



D\ sin (j) sin 29} 



dE 



6(E - £±). (20) 



where are the energies of the bound Andreev states. 
To calculate the current density, information about the 
bound state energies as well as the derivative of the en- 
ergy with respect to phase difference is thus needed. The 
bound state energies are given by the-aeros of the denom- 
inator Z [Eq. (|lq)1 at e = 0, namelyc3 



cos 29 — R cos x + D cos < 



(21) 



The energy of the Andreev levels as a function of phase 
difference cf> is plotted in Fig. ||. In the figure it is shown 
that the Andreev levels appear in pairs, labeled by n, 
with an upper (+) and a lower (— ) level (referring to 
E > 0). The index n is zero for the pair of levels with 
positive energy closest to Ep. In the case of one single 
bound state, the level is labeled by Eq . 




< 



and splits asymmetrically between the two horizontal 
arms 2 and 3, 

2e e 2 

iin.j2,3 = ± — — {1 - COs(2# - fyj). + COS 0(cOS X COS 26)} 



h Z 




(19) 
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From the relations i + (E) = —i + (—E) and i~(E) = 
i~(—E) one can calculate the current densities for all 
energies inside the gap \E\ < A. The continuum current 
density, for energies outside the gap \E\ > A, is calcu- 
lated in the same way. However, since the Andreev reflec- 
tion probability decays very rapidly outside the gap, the 
Andreev resonances become very broad and contribute 
much less to the current. Only the quasiparticles injected 
from the superconductors contribute significantly to the 
current, as will be discussed below. The full formulas for 
the continuum current density for a symmetric junction 
I = is presented in Appendix A. 



<t>/TC 

FIG. 4. Andreev bound state energies as a function of 
phase difference <f) for different lengths L — (left), L ~ £o 
(middle) and L ^> £o (right) of the junction with D — 0.7. 
Solid lines are for a symmetric junction / = 0, dashed for an 
asymmetric one. A gap opens up in the spectrum at <f> — 
due to the asymmetry. 

The derivative of the bound state energy with respect 
to phase is obtained by differentiating Eq. (|2l|), giving 



dm 



D sin0 
2 sin 26 
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A 2 - (E, 



±\2 



L 
Tlvf 



l R sm X 
hv F sin 20 



(22) 



The expression for the sum current density is given by 
inserting Eqs. @-(E|) into Eq. (pj), giving 



2e^ r dE 



eh 



(23) 



where the relation sgn[(<iE /d<fi) sin tfi sin 29] = — 1 [see Eq. 
(B9)] has been taken into account. The expression ( |2^ ) 
coincides with the equation for the Andreev bound state 
currentEj derived directly from the BdG equation. From 
the alternating slopes of the energy-phase relation E(<fr), 
plotted in Fig. [|, it is clear that the sum current density 
(~ dE/d(f>) changes sign between two subsequent An- 
dreev resonances (see Fig. ||). 

The anomalous current density i a is given directly by 
inserting Eq. ( |20| ) into (|l7|), namely 



2e . 
-a— sgn(sm 




6(E - Et). 



(24) 



For a symmetric junction I = 0,cosx = 1, the anoma- 
lous current density does not change sign as a function 
of energy, opposite to the sum current density (see Fig. 
|3|). For finite asymmetry, the anomalous current might 
change sign. However, this does not lead to strong sup- 
pression of the total anomalous current, as will be shown 
below in section VB. 

The injection current ii n j = i^ — is proportional to 
e 2 and thus goes to zero for e<1. We approximate the 
injection current in the weak coupling limit by the first 
order term in e, given by inserting the expression for e/Z 
in the zero coupling limit into Eq. (EOh 



h 

dE± 



n,± 



X + D(cosx — cost 



I sin 201 



S(E-Et). 



(25) 



The injection current density is closely related to the 
anomalous current density i a , in the sense that the injec- 
tion current density is positive for all energies and values 
of the phase difference 4>. 



C. Structure of the nonequilibrium current 

Including the continuum contribution from the super- 
conductors (Appendix A) in Eq. (Q) , we can finally write 
down the structure of the total current in each lead: 



dE 



n h ) + i s n F 



, (26) 



where i s is the current density from the quasiparticles 
injected from the superconductors. The equilibrium cur- 
rent (V — 0) flowing in leads 2 and 3 is given by 



I eq = I dE [i + + i s ] n F 



(27) 



while in lead 1 it is zero. Subtracting the equilibrium 
current from the total current we get the nonequilibrium 
current in the horizontal leads 2 and 3. We divide the 
nonequilibrium current into the the regular current I r 
associated with the nonequilibrium population of the ex- 
isting resonant states, 



dE 



— (n e 
2 y 



n h - 2n F ) 



(28) 



the anomalous current I a associated with the essential 
modification of the Andreev states due to the open nor- 
mal lead, 



dE 



and the injected current I\ 



h 



dE 



(29) 



(30) 



With these definitions, the total currents in leads 2 and 
3 may be written as 



h — Ieq 



Ir ~t~ la Iinj, 2; 
Ir ~t~ la ~1~ Iinj ,3 ; 



(31) 



where I, 



hnj,2+Imj.3- As discussed in Section II, the 
separation of the anomalous current is arbitrary, and has 
physical meaning only in the weak coupling limit when 



0. 



In the weak coupling limit, the integrals in Eqs. (£ 
become sums over resonant states 



I, 



ly. 



dE± 



[n e (Et)+n h (E±)-2n F (Et)], (32) 



-a— sgn(sin < 



E 



cos x — COS ( 



|sin20| 



dEt 



[n*(Et)-n h (E±)]. (33) 



The equilibrium current for energies \E\ < A is given by 
inserting Eq. @ into d27|), 



(34) 



2e \ ^ dE^ , ±s 



n,± 



For energies above the gap, the equilibrium current re- 
sults from quasiparticles injected from the superconduc- 
tors only, since this current is the only continuum current 
being finite in the weak coupling limit (see Appendix A). 
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IV. JOSEPHSON CURRENT OF A SHORT 
JUNCTION 

For a short junction L = I = 0, there is exactly one 
resonance for positive energies < E < A. For no cou- 
pling to the normal reservoir e = 0, this resonant Andreev 
state is converted into a bound Andreev state, with the 



dispersion relation E = Ay 1 — D sin 2 '((f) / '2) . The equi- 
librium cuKccnt of a short junction is thus given by the 
well knownEa relation 



Ieq 



-:A 



D sin 4> 



2 A /l-Dsin 2 (0/2) 



:tanli(JSo /2fcT). (35) 



The continuum current is zero, which can be seen by 
putting L — (/3 = 0) in the equations for the contin- 
uum current in Appendix A. At zero temperature and 
zero applied bias, only the level with negative energy 
— Eq is populated. For an applied a voltage bias V > 0, 
the electron (hole) population is shifted upwards (down- 
wards) in energy. When the voltage eV = Eq , the en- 
ergy of the resonant level, the level becomes populated 
and there is an abrupt jump of the current. The regular 
part of the current, I r , jumps an amount SI r = —I eq , 
thus cancelling the equilibrium Josephson p.uxrent. This 
has recently been observed in experiments. Bui 
The anomalous current jumps by the amount 



_ 0.4 
< 

~ 0.2 



0.5 1 

eV/A 

FIG. 5. The currents I eq + Ir (dash-dotted), I a (dashed) 
and the total Josephson current I eq + I r + I a (solid) in the 
horizontal leads 2 and 3 as a function of voltage V at T = 
for a short junction L — with D = 0.8, <f) = 3tv / 4, e — 0.01 



and a = — 1. 

for eV > E- 



The total current is I eq for eV < E n and SI a 



The effect of finite temperature in a zero length junc- 
tion is merely to smear the steps in the IVC. 

In the symmetric case (I = 0) it is interesting to extend 
the discussion to a longer junction with two resonant lev- 
els (see Fig. ^) , since the cuxpnt distribution between the 
levels becomes nontrivial.OEj In the limit D <C 1, both 



resonances have energies close to the gap edge, E w A, 
and with the additional approximation /3/2 > \fi5 we 
obtain the expression for the derivative of energy with 
respect to phase [see Eq. (p2|)] 



dE* 



± 



A-/D L 



sin(0) 



Co 



(36) 



sin(0/2)U/l -D sin 2 (0/2) 



The equilibrium bound state current becomes propor- 
tional to I h eq ~ <IEq /d(f> + dEq /defy ~ D (taking terms of 
order D into account), but for the currents of the individ- 
ual levels ~ \[T). The resonant levels thus carry opposite 
"giant" currents which almost cancel in equilibrium. For 
L > 0, we also have to take the continuum contribution 
into account. In has been shownLH that the continuum 
contribution to the equilibrium current is I° q — —l/2I^ q> 
thus giving the total equilibrium current I eq = \j2l\ q . 

At zero temperature, when a voltage equal to the low- 
est lying level eV = Eq is applied, the regular and 
anomalous current jumps 



SI r 



eA L 



L>sin(0) 



k & 2|sin(0/2)U/l -13 sin 2 (0/2) 



SI a 



lA L 



h & ,/l-Z3sin 2 (0/2) 



(37) 



(38) 



Both jumps are proportional to yD, and the magnitude 
of the total current at Eq < eV < Eq is then much 
larger than the equilibrium current. 

When the bias voltage is further increased to eV = Eq 
there is a second current jump: the regular current jumps 
in the opposite direction and becomes equal to the small 
negative bound state equilibrium current I r — — 1\ . The 
anomalous current, however, again jumps 8I a in the same 
direction. For voltages eV > Eq the total current in the 
junction is thus I^ q + 281 a . The full formulas for all the 
individual currents including temperature dependence is 
given in Appendix C. 



V. JOSEPHSON CURRENT OF A LONG 
JUNCTION 

We now discuss the Josephson current in a long (L 3> 
Co) symmetric (I = 0) junction, and we treat the effects 
of asymmetry below. In a long junction there are man - 
(N = [L/(£on)]) pairs of resonances, as seen in Fig. 
The width of each resonance is T — ehvF / (2L). This 
width must not be too small if the quasiparticles are to 
be able to enter and leave the junction without being 
scattered inelastically (r > Tivp/li). This gives an upper 
limit for the length L < ^e. 
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The derivative of energy with respect to phase dEldS 
in Eq. (p2h, which determines the current in Eqs. (|32|)- 
(B4h, can be simplified in a long junction L ^> £ , 



±- 



L>sin(0) 



(39) 



4|sm(0/2)U/l-rJ>sin 2 (</>/2) 



This expression holds everywhere except close to the gap 
edge, A — E n ~ (Hvf/L)(^o/L), a distance much smaller 
than the energy distance nhvp/L between the pairs of 
levels. Therefore, equation ( |39| ) can be used for calcula- 
tion of the currents of all levels except the last pair of 
levels closest to the energy gap. The current from this 
last pair of levels must always be treated on a separate 
footing. 

According to Eq. (39), each of the Andreev levels carry 
a current of the order of l/L. Furthermore, as follows 
from the exact Eq. (^H), each pair of levels carries a 
small net current, dE+/d<j> — dE~/dcf), of the order of 
(l/L) 3 . The sum of the currents from all bound states 
is thus determined by the current (~ l/L) from the last 
pair of levels. 



A. Equilibrium current 

For the equilibrium current of a long junction, the con- 
tribution from the Andreev bound states at \E\ < A and 
from the continuum at \E\ > A are of the same magni- 
tude (see Fig. |6|). The continuum current (see Appendix 
A) is given by 



e 
h 



-A 



dEnp 



4D sin <j) sin (3 sinh 2^/ c 



(cos (3 cosh 27 c — R — D cos 4>) 2 + (sin (3 sinh 2^ c 



(40) 



with-^jy c = arccosh(EyA). —Following the method by 
Ishiic£l and Svidzinsky et alE3, one can rewrite this inte- 
gral as a sum over the residues, 



r: 



-I h eq + AkTm^i(E p ), 



(41) 



where the first term results from the poles of the current 
density i(E) in Eq. (|f|) and the second term from the 
poles of the distribution function np(E), given by E p = 
i2kTn(l/2 + p). The first term in Q is the current 
carried by the bound states with negative sign. The total 
equilibrium current I eq is then just given by the second 
term. For zero temperature this sum over poles turns 
into an integral and the total equilibrium current is then 
given by 



I eq (T = 0) = - 



Hvf \D sin(0) arccos(i? + D cos < 



h L 



2tt\ sin(^/2)| v/l - L>sin 2 (0/2) 



(42) 



At high temperatures hvp/L <C kT <§; A, only the the 
first term (p — 0) in the sum ( [4l|) needs to be included, 
and the equilibrium current becomes 



I eq (kT^hv F /L) 



eA /4fcT 



Dsmde- 27rL/iT , 



(43) 



where £t = Tivp/kT. The equilibrium current for a 
long junction, at finite temperature is thus exponen- 
tially Sfflall.EZl. Expressions (^2|) and (^) extend earlier 
resultal5E3 to the case of arbitrary transparency D of the 
junction. 



0.4 



< 0.2 




0.02 




0.01 




■ 




0.01 







kT/A 



0.1 



FIG. 6. The equilibrium bound state (a) and continuum 
(b) currents and their sum (dashed) as a function of length 
L for finite kT = 0.2A, D = 0.8, (f> = 3tt/4 and e = 0.01. 
There is a cusp in both the bound state and continuum cur- 
rents when a new bound state forms out of the continuum. 
The total equilibrium current, however, dies monotonically 
with increased length. Inset: The equilibrium bound state 
and continuum currents as a function of temperature for a 
long junction L = 15£o with D = 0.8, <j> = 3tt/4 and e = 0.01. 
The bound state current (a) decreases from Ie q (T = 0) = i* 
to —leqi when the temperature is increased from zero to 
kT <^ Tivf/L. The continuum current (b) is unaffected in 
this temperature regime. 

We are also interested in analyzing the separate be- 
havior of the bound state current, because this current is 
revealed in nonequilibrium, as will be discussed in detail 
below. Therefore, using relation (39) we can write Eq. 
(M) on the form 



I b = 



e hvi 



D sin((/>) 



'* L 2| wx(<l>/2)\y[l -Dsin 2 (0/2) 

/v-i 

[ta,nh(E- /2kT) - taah(E+ /2kT)] 



n=0 

+i*tanh(A/2fcT). 



(44) 
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The term i* results from the last pair of levels at E w A, 
and is of the order 1/L, as discussed above. At T = 0, 
the sum in Eq. (]44|) is zero, and we thus find that 
Ig q (T = 0). When the temperature is increased, 



the sum in Eq. (44) starts to contribute with nega- 
tive sign and the bound state current is decreased. The 
continuum current (and also i*), however, is indepen- 
dent of temperature for kT <C A, since it is an in- 
tegral over states with \E\ > A [see Eq. (40)]. At 
kT hvp/L the total equilibrium current is exponen- 
tially small (see Eq. and has thus decreased an 

[kT > hv F /L) » I eq (T = 0). 



amount I eq (T = 0) 



This is thus solely due to decrease of the bound state 
current, as shown in the inset in Fig. |[ 



B. Regular current 

The regular current can be written, inserting relation 
(EM) into Eq. (p2), on the form 



e hvi 



Dsm(4>) 



h 2L 



2| sin(<t>/2)\<Jl - D sin 2 (0/2) 



N-l 



Y / [g(E-)-g(E+)] + -g(A) 



(45) 



n=0 



where g(E) = tanh[(E + eV)/2kT] + tanh[(£ - 
eV)/2kT}-2 t&nh(E/2kT). The regular current I r jumps 
up or down every time eV = E^ [see Fig. [7j. Each cur- 
rent jump has the magnitude 



e hv F 
SIr = h — 



D sm{4>) 



2|sin(0/2)|y / l - Z) sin 2 (0/2) 



(46) 



FIG. 7. The equilibrium current I eq plus the regular cur- 
rent vs voltage. L = lO^o, 4> = 7r/2, D = 0.8, e = 0.05. 
Solid line - T = 0, dashed-dotted - fcT = 0.04A, dashed 
- fcT = 0.07A. The regular current jumps alternating by 
±SI r every time the voltage is equal to the energy of an An- 
dreev resonance. For kT ^> Kvf/L and eV > A the current 
T A- T — T c 

leq \ J-r — J-eq 

It is interesting to study the sum I eq + I r , plotted in 
Fig. 0, at temperatures kT fivp/L. In this tempera- 
ture regime the equilibrium current is exponentially small 
and also the regular current steps in the IVC in Fig. 
are suppressed. For a voltage eV ~ A, the last level, 
carrying the major part (i*) of the bound state current, 
is populated and the current I eq + I r jumps to I° q , the 
value of the continuum current, since all bound states are 
populated. This current I£ q is of the order of 1/L and 
the current I r + I eq is increased from zero to be ~ 1/L 
when increasing the voltage from V — to eV ~ A. 



C. Anomalous current 

The anomalous current is given by inserting Eq. 
into Eq. (|3S" 



[h(E+)+h(E-)], (47) 




AL 1 -Dwa?(6) ^ 

v™ ' n—0 



where h{E) = ta,nh[(E-eV)/2kT}-ta,nh[(E+eV)/2kT}. 
We have neglected the current from the last level close 
to E = A, because the currents of all levels add up and 
the current from the last level is negligible. The IVC at 
zero temperature looks like a staircase, as shown in Fig. 



at zero temperature. At voltages eV > A, the regular 
current is the sum of all states in the range < E < A, 
and is equal to the negative bound state current —i*. 



< 



+ 




-1.0 





eV/A 



1.0 




FIG. 8. The anomalous current as a function of voltage V 
for (a) (j) = tt/4 and (b) 4> = 3tt/4 for L = 10£o, £> = 0.8, 
e = 0.05 and a — — 1. Temperature T = (solid) and 
T = 0.1 A (dashed line). The current steps with magnitude 
5I a are smeared to a straight line for kT 3> frvp/L. Up- 
per inset: The critical anomalous current at eV = A as a 
function of transparency D for coupling constant e = 0.1. 
Due to finite coupling e, the critical current always goes 
to zero for R = 0. Lower inset: The anomalous current 
I a (eV = A, kT S> Kvf/L) as a function of phase difference <j> 
for different transparencies D = 0.1, 0.5 and 0.9. The highest 
amplitude corresponds to the highest transparency and vice 



The magnitude of the current step at zero temperature 
is given by 



Sla 



e hv} 



h 2L 1 - D sin 2 (0/2) 



(48) 



At temperatures larger then the interlevel distance, 
kT ^> frvp/L, the staircase IVC is smeared out to a 
straight slope, as shown in Fig. |^. The exact position of 
each level becomes irrelevant and we can write the sum 
over bound states in (ft7| ) as an integral, noting that the 
expression dE/dn — itTlvf/L holds for all levels in the 
sum @, 

N 

^2[h(E+) + h(E-)} 

n=0 

or /"A 

w / d£[tanh(£ + eV) - tanh(£J - eV)] 

nhv F J 



4L 



f(V,T) 



(49) 



where 



and the anomalous current takes the simple form 



MT). 



hir[l -L> sin 2 (0/2)] 



(51) 



In the limit fiv F /L < kT < A, f(V,T) = mm(eV,A): 
the anomalous current thus scales line arly with applied 
voltage up to A. It follows from Eq. ( pl|) that / a is in- 
dependent of the length of the junction, being the sum 
of N ~ L levels which each carries a current /„ ~ 1/L. 
This gives that the anomalous current roughly is equal 
to the total equilibrium current of the short junction. 
The critical anomalous current is plotted with respect to 
transparency in the inset in Fig. ||. In the limit D -C 1 
it is given by (J ) c = (e/h)(y/D/ir)f(V,T). It is propor- 
tional to the first power of A for T close to T c , therefore 
surviving up to kT w A. The anomalous current-phase 
relation (see inset in Fig. ||) is 27T periodic and resembles 



that of the equilibrium Josephson current. The direction 
of the anomalous current is however proportional to a, 
i.e dependent on the phase of the scatterer at the con- 
nection point, which is not the case for the equilibrium 
Josephson current. 

To get the complete picture of the Josephson current in 
a long junction, / = I eq + + I a is plotted as a function 
of bias voltage for different temperatures in Fig. pi 
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eV/A 



FIG. 9. The total current I = I eq + I r + I a as a function of 
voltage. At zero temperature we have + I eq (dash-dotted), 
I a (dashed) and the total current I r + Icq + la (solid) . The 
total current for temperatures kT S> Tivf/L is plotted (dot- 
ted). Junction parameters are D = 0.8,0 = 37r/4, L = 20£o, 
e = 0.01 and a = -1 



The zero temperature total Josephson current oscil- 
lates strongly around a constant slope as a function of 
voltage, showing steps whenever the voltage passes an 
Andreev level. The step structures are washed out for 
temperatures kT 3> Iivf/L, and in this limit the total 
current roughly coincides with the anomalous current, 
given by Eq. ( |5l| ) . 



D. Asymmetric junction 

The effect of asymmetry is most pronounced in the 
long limit when the asymmetry is much larger than the 
coherence length but much smaller than the total length 
of the junction, L I I n this limit, the deriva- 

tive of energy with respect to phase dE/d(f> in Eq. |2^ 
reduces to the expression of a symmetric long junction 
@, since |sin20| > -R|sinx| (see Appendix B). The 
equilibrium current I eq and the regular current I r are 
not substantially changed in comparison to the symmet- 
ric case. In contrast, the anomalous current is modified 
in a nontrivial way, taking the form 

n L 
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N 



St 



cos x — cos ( 



n,± 



— (D cos 4> + R cos x) 2 



(n e ~ n h ), (52) 



obtained by inserting Eq. @ into @. For T = the 
step structure in the IVC is modified due to the change of 
Andreev levels as a result of the asymmetry (see Fig. ^) . 
Already for small asymmetry I ~ £o the anomalous cur- 
rent might change dramatically. Depending on the phase 
difference of the junction, the IVC is renormalized and 
changes sign for — tt/2 < <f> < ir/2. 



< 




eV/A 

FIG. 10. The asymmetric anomalous current I a vs voltage 
for different asymmetries (a) I — 0£o, (b) I = 2£o and (c) 
I = 40£ for kT > hv F /L, D = 0.8, e = 0.05, L = 60£ , 
<f> — 7r/4 and a = — 1. The IVC is changed dramatically 
already for as small asymmetry Z ~ , if the phase difference 
-tt/2 < 4> < 7r/2 

When the temperature is increased beyond the inter- 
level distance kT ^> hvp/L, the step structure becomes 
smeared and we get a periodic modulation of the IVC on 
the scale of eV ~ %vp/l. This modulation arises from 
the factor cosx- 

When the temperature is further increased to kT ^> 
Tlvf/1 this periodic structure is smeared out and the IVC 
once again becomes a straight line, but with renormalized 
slope. In this high temperature limit the amplitude of 
the terms in the sum in Eq. ( |52| ) oscillates with a period 
Kvf/1- During this period, the filling factors n can be 
taken to be constant, and we can sum over one period to 
get the average value. Performing this summation in the 
continuum limit, we get 



£ r 



cos x — cos ( 



one period 



TlVj- 



(D cos 4> + R cos x) 2 
cos(x) - cos(0) 



over the filling factors following the procedure from the 
symmetric case (49) 



E 



(n e - n h ) 



4/ 



TlVp-K 



f(V,T). 



(54) 



averaged periods 

The anomalous current then becomes 



In 



e D . 
-cr2— ■ — = sm< 
^ TTyR 

|sin(0/2)| 



|cos(</>/2)| 



yi-£cos 2 (</)/2) 



f(V,T), 



(55) 



which is independent of both the length L and the asym- 
metry I, We also find that the renormalized anomalous 
current becomes 7r-periodic. This can qualitatively be 
explained by the fact that the 27r-periodic part of the 
anomalous current density is very sensitive to asymme- 
try, oscillating fast with energy on the scale of Hvf/1, be- 
coming washed out during summation over bound states 
at high temperatures kT ^> Uvp/l- The it periodic part 
of the current does not have this sensitivity and is the 
only part of the anomalous current that survives. The 
asymmetric anomalous current-phase relation is shown 
in Fig. H 




FIG. 11. The asymmetric anomalous current I a at eV = A 
and kT 3> Tivf/L as a function of phase difference <f> for dif- 
ferent transparencies D = 0.1, 0.5 and 0.9. Inset: The critical 
anomalous current (I a )c, for eV = A and kT ^> Hvf/L, as a 
function of transparency D for coupling constant e = 0.1. 



21 J 1 - (D cos 4> + R cos%) 2 
I I / |sin(0/2)| 



|cos(eV2)| 



The 7r-periodicity and the zeros at (f> = nn/2 give the 
condition that the slope of the IVC must change sign due 
to asymmetry in the range —ir/2 < 4> < tt/2, as shown in 
Fig. [l^. The critical asymmetric anomalous current as 
a function of transparency D is shown in the inset in in 
Fig. [ll]. The critical asymmetric anomalous current as a 
function of transparency D is shown in the inset in in Fig. 
This quantity is energy independent and we can then sum [ll]. The behavior is very similar to the critical anomalous 



I 8RVD \ ^l-ZWO/2) /i_ ,0^(0/2) 



(53) 
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current in the symmetric case, the main difference being 
that the amplitude is reduced by roughly a factor of two. 



VI. INTERFACE BARRIERS 

In any realistic experimental situation, normal reflec- 
tion at the NS-interface, modeled by a barrier with re- 
flection amplitude rb, must be taken into account The 
general expression, considering both the interface barri- 
ers and the midpoint scatterer, becomes analytically in- 
tractable. We can however analyze the case where the 
midpoint scatterer is absent (R — 0) to get an under- 
standing of the effect of NS-barriers on the junction prop- 
erties, and then treat the general case with injection and 
midpoint scatterer numerically. 

In the absence of the superconducting leads (a NININ- 
junction), the two barriers give rise to normal Brcit- 
Wigner resonances for the electrons and holes. Under- 
stating the properties of these resonances turns out to be 
crucial for describing the behavior of Andreev levels and 
current transport. The energies of the electron and hole 
resonances are calculated straightforwardly 



E'„ 



E m = 2Ep 



1 - 
1 - 



7r(n — v) 



kpL 
ir(m — v) 



kpL 



(56) 



where r& = ^/Ri,e lU7T and n{m) are integers denoting the 
index of the electron (hole) resonances. The intersection 
between electron and hole resonances (E^ = E^J is given 



by L 



n+m 



\p / 4(m + n — 2i>) with the Fermi wavelength 



\p = 2n/kp. 

Fig. m 



These normal resonances arc plotted in 



Ln+m-1 



L n+m- 1 
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t L n + m Ln+mf 
FIG. 12. The Andreev levels (solid) and the normal elec- 
tron and hole resonances (dotted) as a function of length L of 
the junction with four Andreev levels in (a) weak resonance 
limit Rb -C 1 (b) strong resonance limit Rb ~ 1. The lengths 
of two subsequent intersections of normal resonances L n +m 
and L n + m -i are shown with arrows 



For the junction with superconducting leads, one can 
in the same way as before calculate the equationJpr the 
bound Andreev states (\E\ < A), with the resulted 

Dl cos c/> + 2R b cos (3 - cos(2 7 - 0) - R.% cos(2 7 + 0) 
-4i? h sin 2 7 cos(A)) = (57) 

where we have defined (3o = ±2E e,h /(Tlvf / L) and +(— ) 
denotes hole(electron) resonance energies. One can draw 
some qualitative conclusions on how the Andreev levels 
are related to the normal resonances by looking at Fig. 
|l2j In the limit of high barrier transparency Rb -C 1, 
the Andreev levels are weakly modified by the barriers. 
In the opposite limit Rb ~ 1, the Andreev levels get 
pinned at the normal resonances, but there are no level 
crossings at the points where the normal electron and 
hole resonances intersect. 

We find the same interlevel distance hvpn/L in the 
junction with the superconducting leads (SINIS junction) 
and normal leads (NININ-j unction). The main difference 
is that the normal resonance move very quickly through 
the junction when the length L increases, while the An- 
dreev levels oscillate up and down. 

Considering Andreev state energies close to the Fermi 
level, E—.<^ A, one can derive a simplified dispersion 
relational 



sin 2 (/?/2) 



J 2 cos 2 (0/2) + 4fl fc sin 2 (/W2) 
(l + i? 6 ) 2 



(58) 



Using this relation we can study the bound state current 
in different length limits. 

In the short limit, L <C £o there are two cases to be con- 
sidered. For nearly transparent barriers Db ~ 1, and thus 
broad resonances T = D^ftvp/L 3> A, one can neglect 
dephasing (putting (3 = 0) and just get the total trans- 
parency of the junction D = D 2 b l{D 2 b + 4R b sin 2 (/3 /2)) 
to be put into the standard zero length junction equilib- 
rium current formula. In the strong barrier case D b <C 1 
the resonances are sharp r <C A and one can not neglect 
the dephasing. Assuming that the resonance is close to 
Fermi !e^eijgy E e ' h <C A, we can put (3 <C 1 in (58) and 
obtai 



E = ± 



yl^ 2 cos 2 (0/2) + {E^ h f 



(59) 



When the resonance is exactly at the Fermi energy 
E e ' h = 0, the Josephson current is given by 



er 

— sin((/>/2) tanh 



rcos(0/2)\ 
2kf J 



(60) 



The critical current at low temperatures (kT <C T) is 
thus smaller than the critical current of a short, clean 
junction by a factor T/A. 

For a long junction L»(o we can calculate the deriva- 
tive of energy with respect to phase, 
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dE 



= ±- 



hvi 



D b sin < 



2L ^/(l + R b )±-[D b cos ( 



AR b cos(/3 )] 2 

(61) 



In the weak barrier limit ft « 1, this just causes oscil- 
lations with length around the clean junction (R b = 0) 
result. In the strong barrier limit Rt, ~ 1, one can dis- 
tinguish two limits: When the length of the junction is 
far away from the length L n+m where the electron and 
hole resonances intersect, the junction is out of resonance. 
The second term in Eq. ( |6l| ) is negligible and the current 
from the individual levels thus becomes 



±— D 6 sm< 



(62) 



In 



It is proportional to D b and thus strongly suppressed, 
the opposite limit, when the length of the junction L = 
L n +m = Xf (m + n — 2v) /4, the junction is in resonance. 
When n + to is even we get the current carried by each 
level 



I = ± 



evF D b smi 



and when 



L 4|cos(0/2)| 
to is odd we get 



± 



evp D b smcj) 
L 4| sin(0/2)| ' 



(63) 



(64) 



We see that the current is proportional to D bl just as ex- 
pected for the junction in resonance. The current carried 
is thus of the order of the single barrier junction current. 
An interesting feature is that the current is dependent 
on the parity of the sum of the electron and hole reso- 
nance indices n + m. When the third lead is connected to 
the junction, the scattering at the connection point just 
splits the Breit Wigner resonances, and the qualitative 
picture for the bound states derived without the third 
lead connected survives. 

To calculate the total equilibrium, regular or anoma- 
lous current, the currents carried by all individual lev- 
els have to be summed up. In the weak barrier limit 
R b -C 1 we just find that all properties calculated above 
for the symmetric junction without barriers hold, with a 
small length dependent modulation ~ R b with a period 
SL ~ Xp. In the strong barrier limit R b ~ 1, the re- 
sult will depend on whether the junction is in or out of 
resonance. 



H 2 




L/A F 

FIG. 13. Short segment of the equilibrium bound state cur- 
rent as a function of length L, to illustrate the resonant be- 
havior. The junction is long L 2> £o with Rt = 0.9, e = 0.01 
and <j> — 7r/2. 

Fig. [l^ shows the resonant behavior of the equilib- 
rium current as a function of length. The current has a 
peak around lengths L = Aj?/4(to + n + 2v). The phase 
dependence of the current at the resonant peaks is well 
described by the expressions for the single level currents 
©and®. 

The anomalous current is also strongly length depen- 
dent and when the junction is in resonance we have an 
anomalous current I a ~ <rD b v 7 RD while when we are out 
of resonance I a ~ aD b \/RD. It turns out that there is 
an anomalous current even without scattering at the con- 
nection point, but it oscillates around zero as a function 
of length with the period ~ Xf- 



VII. INJECTION CURRENT AND 
CONDUCTANCE 

Although the nonequilibrium Joscphson current is at 
the focus of this paper, the injection current that flows 
between the normal reservoir and the SNS-junction is 
also of great interest: it determines the conductance of 
the circuit. The conductance of pNS-structures has been 
studied intensely in recent yearalj and is an interesting 
quantity in itself. It can also be used to determine the 
direction of the Josephson current in the junctior£3 or to 
detect a large Josephson current in the superconducting 
loop which changes the applied external flux vs phase 
dependence, thus modifying the phase dependence of the 
conductance.Ej 



A. Injection current 



We start by discussing the symmetric junction I = 0, 
and comment on the modifications due to asymmetry 
below. The current injected in lead 1 for energies \E\ < A 



13 



can be calculated to lowest order in e by inserting Eq. 
© into Eq. © 



h 



I oos(0/2)| 



v / I?yl- J Dsin 2 (0/2) 



E 

n,± 



[n e (E±)-n h (E±)] 



(65) 



This current is proportional to e (unlike the Josephson 
current discussed above), which is also true for the con- 
tinuum contribution. For \E\ > A, the injection current 



density in Eq. (|3C|), ii n j — i\ = — is obtained 
from Eq. ( |A3| ) . This current density is roughly described 
by the normal current density Aee/h, with oscillations 
around this value due to the resonances in the scatter- 
ing states (see Fig. [l|). These oscillations are strongest 
around E ~ A and decrease with increasing energy. 

The injection current is proportional to the modulus 
\dE/d<j)\, just like the anomalous current, as discussed 
above. The IVC thus has the shape of a staircase, as 
shown in Fig. |lj. 
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FIG. 14. The injection current in lead 1 as a function of 
voltage for (a) <j> = 7r/4, (b) <f> — 37r/4 and (c) <f> — it. Zero 
T (solid lines) and T = 0.05A (dashed line) with D = 0.8, 
L = 10£o and e = 0.05. For eV > A the IVC approaches the 
value of a normal junction 

We can derive expressions for the injection current in 
different length limits. A complete set of formulas are 
given in Appendix A. Here we only present the heights 
of the current steps in the IVC:s in some representative 
cases. 

In the limit of zero length of the junction [L = 0) 
there is only one Andreev level (for < E < A) and the 
magnitude of the step is 



_ e VP |sin(^/2)|cos 2 (^/2) 
1 ~~ e h 2 1 -L>sin 2 (0/2) 



(66) 



Mi = e- 



eA L 



cos 2 (0/2) 



h ^ x /l-rJsin 2 (0/2) 



(67) 



In the long junction limitfL 3> £o) the current is given 
by putting Eq. (§|) into @, 



h =e 



e hv F cos 2 (0/2) 



N 



Y^[h(E-) + h(E+)], (68) 



h 2L l-^sin 2 (0/2)^- o 
The current step at zero temperature is 

^ _ e Hvf cos 2 (0/2) 
1 _ € h L 1 -Dsin 2 (0/2)' 



For high temperatures kT 3> Hvf/L, the sum (|68| ) can 
be converted to an integral, just as for the anomalous 
current (fl9]), which gives the current for eV < A 



(69) 



h =e 



4e cos 2 (0/2) 



h 1 - Di 



f(V,T). 



(70) 



where f(V,T) is given by Eq. @. The jyC thus be- 
comes linear for eV < A, with the slope independent of 
length and temperature, as is seen in Fig. [l4| All infor- 
mation about individual Andreev levels is washed out. 

The effect of asymmetry on the injection current is 
drastic in the limit of a long junction with large asym- 
metry L I £o, just as for the asymmetric anomalous 
current. This shows the strong relationship between the 
two currents. The injection current is given by inserting 
Eq. (H) into (|3C 



% N 

h = € X 



"h L 



n=0 



D sin 2 (0) + R sin 2 x 
(Dcos0 + i?cosx) 



n h ). 
(71) 



Averaging over periods and summing up the filling fac- 
tors, just as in the case of the anomalous current [see Eq. 
([53|)], the injection current becomes 



h 



|sin(0/2)| 3 



|cos(0/2)| £ 



Vl-£>cos 2 (0/2) v /i_D 8 in 2 (0/2) / 



f(V,T). 
(72) 



The injection current in this limit does not depend on 
either the length L or the asymmetry I. It is 7r-periodic, 
7i (0+7r) = Ii (0), for the same reasons as discussed above 
for the asymmetric anomalous current. The implications 
of the 7T periodicity for the conductance is discussed be- 
low. 



while for two levels close to A we get (L ~ £q) 
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B. Conductance 



For the conductance, we discuss the whole range of 
coupling parameters e, not only the weak coupling limit. 
The conductance is defined 



/OO • — 
± [cosh- 



2 (E + eV) 
2kT 



cosh 



_ 2 (E-eV) 
2kT 1 



(73) 



At zero temperature T — the conductance for a sym- 
metric junction I — can be writtenO for eV < A, noting 
that ii(E) = ii(-E), 



G(V,<j>) 



Ae 2 



4 cos 2 (0/2) sin 2 



h [(1 - e)cos2<9 - R- Dcos. 



■] 2 + e 2 sin 2 26>' 
(74) 



At a voltage fulfilling the relation 1 — cos 29 = D(l — 
cos 0)/ VI — 2e, which is exactly at an Andreev reso- 
nance, and a phase difference = mod 2tt, the con- 
ductance is G — 4e 2 /h, which is the conductance for a 
perfect NS-interface. This holds for any transparency D, 
length L and coupling constant e. 

For L — we get cos 6* = eV/A to be inserted into 
(H). This gives rise to a peak in the conductance at 
the voltage eV « Eq , the energy of the single Andreev 
resonance. In the long limit (L ^ £o) for E <C A we 
get 8 = tt/2 — eV/E which results in cemductance os- 
cillations as a function of applied voltaget3c3~E3, with a 
period itHvf/L, the distance between the pairs of An- 
dreev levels. This is made clear by takin g th e derivative 
of current with respect to voltage in Fig. [lj. 

The conductance vs phase relation G(V, <f) is also al- 
tered when voltage is applied. For voltages around 
eV/(Tivp/L) « 7i7T, i.e at an energy between the pairs 
of Andreev resonances, the conductance has a maximum 
around (f> « tt and a local minimum at <f> — [apart 
from the absolute minimum at <f> = tt due to the factor 
cos 2 (0/2) in Eq. @]. For eV/(hv F /L) « tt/2 + nir, i.e 
at an energy between the two Andreev resonances in the 
pair, the maximum shifts to <p = and the minimum to 
(f> = 7r [see Fig. ( 15 )] . This behavior has recently been ob- 
served in both ballistictj and quasiballisticEj junction*. 



It has also been predicted for diffusive junctions.E3'Lll 
This voltage dependence of the conductance holds for 
kT < Hvf/L. 




FIG. 15. The conductance as a function of cj> for different 
voltages (a) eV = 0, (b) eV = 0.075A and c) eV = 0.15A at 
temperature kT = 0.01A. D = 0.9, e = 0.05 and L = 20{ - 

For the additional condition zero voltage V = 0, the 
expression ( [74]) reduces to 



G(0,4>) 



,2e 2 



2cos 2 (0/2) 



h i? + Dcos 2 (0/2)' 



(75) 



For D < R the maximum conductance has a universal 
magnitude G max (Q) — Gne/(1 — e) 2 where Gn = e4e 2 //i 
is the conductance of the junction in the normal state. 

When the coupling is weak, e <^ 1, the maximum con- 
ductance is G m ax ~ cGat, i.e much smaller than the 
normal conductance. In this limit the Andreev reso- 
nances are sharp and there are no available Andreev 
states at Ef, because the scattering at the three lead 
connection point opens up a gap in the Andreev spec- 
trum at <f> = ix (see Fig. ^). The conductance is thus 
strongly suppressed. The conductance rapidly increases 
with voltage and temperature and has a maximum at 
eV, kT ~ Tivp/L. This happens because electrons (holes) 
with energies E > (E < 0) tunnel into the first res- 
onant Andreev state at finite energy. This gives rise to 
a finite energy conductance pea kEl (see Fig. [u]). The 
maximum conductance G max is plotted as a function of 
temperature in Fig. [16] (note that the minimum conduc- 
tance always is zero). 
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FIG. 16. The maximum conductance Gmax/GN as a func- 
tion of T for zero voltage. L = 10£o, D = 0.9 and e = 0.05. 
At T = the maximum conductance is G ma x >C Gn- 
It increases with increasing T, reaches a maximum around 
kT ~ Tivf/L, drops again for fcT > Tivf/L but saturates 
at a constant value G ma :r = Gn for fcT 3> TlVf/L. Inset: 
The conductance as a function of (f> for different tempera- 
tures kT = 0,0.01A,0.025A,0.05A and 0.1A at zero voltage 
V = 0. Increasing temperature from bottom to top at (f> — 0. 
D = 0.9, e = 0.05 and L = 20£ - 

The conductance as a function of phase difference </> be- 
haves differently for different temperatures T, as is shown 
in the inset in Fig. |l6|. For kT 3> hvp/L, a limit only 
accessible for the long junction i > (o, the maxima of 
the conductance around (f> sa n are shifted to a maximum 
at <fi — mod 2ir. This holds independent of voltage ap- 
plied. A similar effect in a quasiballistic system has been 
reported by Dimoulas et alEj 

In the weak coupling limit e « 1 we can use expression 
( f70| ) to get the conductance in the long limit for Tivf /L <C 
kT < A and eV < A 



G(cf>) 



4e 2 



cos 2 (</>/2) 



ft l-L>sin 2 0/2) 



(76) 



which is independent of voltage, temperature and 
length of the junction and has a maximum G max = Gn 
at <fi = mod 2ir. 

In the same limit we get the conductance in the asym- 
metric junction from expression (72) 



i - Vd 



|sin(0/2)p 



|cos(0/2)| E 



v/l-^cos 2 ^) x /l-Csin 2 (</)/2) 



(77) 



It is 7r-periodic and has a maximum for <f> — ir(n + 1/2) 
and a minimum for (f> = nn. This 7r-periodicity can be 
qualitatively explained by considering the lowest order 
paths giving rise to the conductance. 




e w 

FIG. 17. The paths in the asymmetric junction giving the 
first order terms of the conductance. Electrons are drawn with 
solid lines, holes with dashed. The grey ellipse denotes the 
effective scatterer due to the three lead connection. The upper 
paths give rise to a 2tt periodic component of the current, 
suppressed at finite temperature kT 3> Tivf/L. The lower 
paths, time reversed, give rise to a i periodic component of 
the current, not suppressed by temperature. 

The upper paths in Fig. [l7], corresponding to an in- 
jected electron and giving rise to an outgoing hole, pro- 
duce a current of the order ii ~ \ e i ^ L+ ^ + e i ^ R+ ^ | 2 = 
2 + 2cos(0 + %). This part of the current is 27r periodic 
in the phase, rapidly oscillating in energy with a period 
Hvf/1- It is washed out when summing up the levels in a 
long junction at finite temperature. The lower paths in 
Fig. nTl corresponding to an injected electron giving rise 
to an outgoing electron, produce a current of the order 
t! - \d*e l ^ + ^ + d*e l (-*+' 3 )| 2 = D(2 + 2 cos 20). This 
part of the current is it periodic in phase and not sensitive 
to asymmetry. 

The discussion about the periodicity of the conduc- 
tance oscillations with respect to phase goes back to the 
early eighties. A 7r-periodic contribution to the weak lo- 
calization correction to the conductance in a similar sys- 
tem was predicted by Scivak et al.Ej and discussed fur- 
ther by Altshuler et alo It has been shown in numerical 
simulations for a structure similar to ours that the full 
conductance, i.e not only the weak localization contribu- 
tion, might become 7r-periodic at finite temperatures iiil A 
large 7r-periodic conductance oscillation with phase was 
also observed in diffusive samplesE 2 ]. Whether the ex- 
planation to the crossover from 27T to 7r periodicity with 
increased temperature discussed above can account for 
these observations remains to be investigated. 



VIII. FOUR TERMINAL JUNCTION 

One problem with the three terminal junction is, 
as discussed above, that it is not possible to separate 
the injection current from the Josephson current in a 
clear wot ..for arbitrary coupling e. In a four terminal 
junctionE^rEjH'OO this is possible under certain condi- 
tions, which makes it interesting to discuss this configu- 
ration separately. 

We consider two different types of junction configura- 
tions (see Fig. n8|). The upper junction is a straight- 



16 



forward extension of the three terminal device pictured 
in Fig. [l]. Two normal reservoirs are connected to the 
normal part of the junction. The reservoirs are then con- 
nected to the grounded superconducting loop via voltage 
sources biased at V\ and V4 respectively. 
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FIG. 18. Two different setups of the four terminal junction. 
In (a) the normal reservoirs are biased independently at V\ 
and Va with respect to the superconducting loop (grounded), 
in (b) only the potential difference between the normal reser- 
voirs, V, is determined. In the right figure, a close up of 
the junction area is shown, with the direction of the currents 
showed with arrows 

The current injected from a normal reservoir is split 
at the connection point. One part of the current flows 
through the junction directly into the other reservoir and 
the other part of the current is divided between the leads 
2 and 3. In the general case, the currents in lead 1 and 
4 are not equal I\ 7^/4. However, by adjusting the po- 
tentials Vi and V4, the currents in the vertical lead can 
be put equal and from current conservation at the con- 
nection point it follows I2 = I3 and we have a clear 
separation between the injection current I\ = I4 and 
the Josephson current circulating in the superconduct- 
ing loop I 2 = h ■ 

In the lower junction, this separation follows directly 
from current conservation at the connection point, since 
I2 = In this junction a bias V is applied between 
the normal reservoirs, which are only connected to the 
superconducting loop via the four lead connection point. 
We define the potential of the superconducting loop in 
this junction to be zero and the potentials of the normal 
reservoirs to V\ and V4 resepectively, just as for the upper 
junction. Our biasing arrangement then gives V — V\ — 
V4. The condition of current conservation, Ii(Vl, V4) = 
h(Vi, V4,), gives a second condition on V\ and V4. With 
these definitions we can calculate the current in both 
junctions in the same way. 

The cross-shaped connection point is modelled by the 
scattering matrix 



S : 



r± \T e V* rf J 

y/e r d y/e 

\fe d r 

d± v 7 ^ v 7 ^ r± 



(78) 



where the e describes the coupling of the SNS junction to 
the vertical normal lead (0 < e < 0.25). The horisontal 
scattering amplitudes now obey the relations Re(rd*) — 
—e and D + R = 1 — 2e. The same holds for the vertical 
scattering amplitudes r± and d±. 



The current densities 



,.e(h),l(4) 



with the upper index 1 



or 4 denoting the lead from which the quasiparticles are 
injected, are calculated in the same way as in the case of 
the three terminal junction. Due to the symmetry of the 
scattering matrix, quasiparticles injected from leads 1 or 
4 give rise to the same current density in leads 2 and 3, 

■ e(h),l .e(h)A 
i e l 2(3) ~ Z 2(3) ■ 

The expressions for the sum and anomalous current 
densities become very similar to the three terminal ex- 
pressions [see Eq. ( |l6| ) and (p"7|)], i.e one just changes 
e — * 2e [also changing Z = Z(e — > 2e)] and divides by 



two, noting that i + ' = 



and iz = 



Neither the 



vertical transparency D±_ nor the reflectivity R± thus 
appear explicitly in these expressions. The factor one 
half simply reflects that there are two normal leads con- 
nected to the normal part of the junction. In the limit 
of weak coupling e <C 1, the sum of the current densi- 
ties from both normal reservoirs is equal to the current 
density from the single normal reservoir in the three ter- 



minal junction, i 



+.1 



+ i 



+A 



and 



i a (simply 



reflecting that one cannot create more Josephson current 
by adding more normal leads). 

In this weak coupling limit the current in the horisontal 
lead I = I2 = I3 is given by 



1 = 1, 



Ir(V 4 )} 



W)] (79) 



with I eqj I r and I a the same as in the three terminal case, 



given by Eq. (|32j)-(|34j). Noting the relations I r (-V) = 
I r (V) and I a {— V) = —I a (V), we see that (i) for bias 
V\ + V4 = the anomalous current is zero, and (ii) for 
Vi — V4 = the regular current is zero. We can thus 
control the regular and anomalous currents in the upper 
junction in Fig. [l8] independently by adjusting either the 
potential difference V\—Va or the sum V1+V4 between the 
normal reservoirs, keeping the other quantity constant. 



A. Injection current and conductance 

The injection currents I\ and J4 in the four terminal 
device is qualitatively different from the injection cur- 
rent in the three terminal device, since in the four termi- 
nal junction the injected quasiparticles from one normal 
reservoir can travel directly through the junction to the 
other normal reservoir. 

The current in leads j — 1,4 can be written [see Eq. 
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dE 



--,1 



- oo 

• + ,4 —,4 

-^-(r^ 4 W M ) + ^-K' 4 



(80) 



The symmetry of the junction gives that i 4 ' = —i 1 ' 
and ' 4 = — ii' ■ This leads to that for V± = —V4 = V/2 
the currents in the vertical leads are equal, i.e Jj n j = I\ = 
J4 and thus no injection current flows into the supercon- 
ductors. In this case Eq. (80) reduces to 



dE 



-{ri 



e,l 



(81) 



Considering the simplest case with a symmetric junc- 
tion / = without barriers at the NS-interfaces. The 
injection current density is given by straightforward cal- 
culation 



-• 4 = 4 



D 



jr[([l 
Z f 



-R) [R± -D±- 2Re[d(d± - r±)] sin^^/2)] 
-l + 2e + (Dcoscj) + R) cos26>) (R± - D±)]} 



2e] cos 26 — D cos <f> 

sin 2 26> 
(82) 



Some general comments can be made about the conduc- 
tance G = dIi n j/dV. When the vertical and horisontal 
leads are decoupled (e — > 0), the conductance reduces 
to G = (e 2 /h)D±, the conductance of the normal ver- 
tical channel. For finite coupling, an additional term is 
added to the conductance SG ~ e. This additional term 
8G is dependent on the phase difference </>, but it is also, 
unlike for the three terminal junction, dependent on the 
scattering amplitudes r,d,r± and d±, i.e not only the 
scattering probabilities R, D, R± and D±. This becomes 
clear when we note that we can rewrite the expression 
Re[d(d±_ - r±)] = 1/2[R± -D±+ ({R± - D ± ][R-D] - 
4Im[r±d* L ]lm[rd*]) / (1 — 4e)], i.e dependent on (t<j±, just 
like the anomalous current. This contribution can be ex- 
plained qualitatively by interference between quasipar- 
ticle paths where one path describes scattering in the 
vertical direction and the other one in the horisontal di- 
rection, thereafter leaving the junction. 

In the case of zero temperature and voltage, the con- 
ductance becomes 



G 



2e 2 



[Da 



D x - R±+ sin 2 (^/2)Re[d(dj_ 



R + L> cos 2 (0/2) 



(83) 



which is independent of the length L of the junction. The 
conductance at zero phase difference given by G((f> = 0) = 
(2e 2 /h)[D±/(l — 2e)]. From this value, the conductance 
then increases or decreases, depending on the phases of 
the scattering amplitudes, monotonically with <j) — > n, as 
is seen in Fig. na 




FIG. 19. The zero voltage, zero temperature four terminal 
conductance G as a function of phase difference <fi for different 
horisontal transparency D = 0.2 (dashed), D = 0.5 (solid) 
and D — 0.8 (dotted). The product aa± — +1 for increase 
of conductance due to finite phase difference, aa± = — 1 for 
decrease. The vertical transparency is D± — 0.5 and e = 0.05. 

For finite energies, the injected quasiparticles tunnel 
into the Andreev resonances and there is a finite bias 
anomaly of the conductance. However, there is not al- 
ways a peak in the conductance. This is made clear by 
considering the conductance at finite voltage, zero tem- 
perature and zero phase difference (to avoid dependence 
on <to"_l), given by 



1 



{R ± - D ± )[cos 2 6- (1 - 2e) 



(1 - 2e) 2 sin^6» + 4e 2 cos 2 



(84) 



In the weak coupling limit, e <C 1, for R± > D±, 
there is an increase in the conductance for finite ener- 
gies, just as for the three terminal device. For R± < D±, 
i.e for highly transmissive junctions in the vertical di- 
rection, however, there is a decrease in conductance for 
finite energies. The presence of Andreev resonances thus 
decreases the probability of a_quasiparticle to be trans- 
mitted through the junction.E3 

Other properties of the conductance are similar to the 
three terminal junction, also taking the scatting phases 
into account (via a and a±). The conductance is periodic 
in voltage with a period ttJivf/L, the distance between 
the pairs of Andreev resonances, for temperatures lower 
than this inter pair distance kT Tivf/L. For such low 
temperatures the phase dependence of the conductance 
is also dependent on the voltage. For high temperatures 
kT 3> Kvf/L, all the features of the individual Andreev 
levels are washed out and the conductance becomes inde- 
pendent on junction length and applied voltage eV < A. 
For eV > A, the conductance is equal to the normal 
conductance of the junction Gn = (2e 2 /h)(D± + 2e). 
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IX. CONCLUSIONS 

We have analyzed the equilibrium and noncquilibrium 
Josephson currents and conductance in a ballistic, multi- 
terminal, single mode SNS junction. The noncquilibrium 
is created by means of quasiparticle injection from a nor- 
mal reservoir connected to the normal part of the junc- 
tion. By applying a voltage V to the normal reservoir, up 
to the superconducting gap A, the equilibrium current of 
a short junction L <C £o can be suppressed. When the 
junction is longer L > £q, the direction of the Joseph- 
son current changes sign as a function of applied voltage. 
For a junction longer then the thermal length L 3> £t, 
the equilibrium Josephson current is exponentially small. 
The nonequilibrium Josephson current in this regime is 
dominated by the anomalous current, arising from the 
modification of the current carrying Andreev states due 
to coupling to the normal reservoir. This anomalous cur- 
rent scales linearly with applied voltage and saturates 
at a magnitude of the order of the equilibrium current 
carried by a short junction, / ~ eA/h. 

The conductance oscillates as a function of the phase 
difference cf> between the superconductors, with a period 
of 2n in a symmetric junction. The position of the con- 
ductance minima, <fi — or ir, is dependent both on ap- 
plied voltage and temperature. The conductance exhibits 
a finite bias anomaly, at eV ~ hvp/L, the position of the 
first current carrying Andreev level. 

Asymmetric injection gives rise to oscillations of all 
currents on the scale of eV ~ Kvf/1 where I is the length 
difference between the two leads. At temperatures above 
this energy, these oscillations are smeared and a we get 
renormalizcd anomalous and injection currents that are 
7r-periodic. 

Introduction of barriers at the NS-interfaces give a 
strong length dependence of all currents, governed by 
the Breit Wigner resonances between the normal barri- 
ers. There are resonant current peaks at lengths where 
the normal electron and hole resonances cross. 

Connecting a second normal reservoir to the normal 
part of the junction allows a clear separation between 
the injection current, flowing between the two normal 
reservoirs, and the Josephson current, flowing between 
the superconductors. 
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APPENDIX A: CONTINUUM STATE CURRENT 

Here we present formulas for the continuum current for 
a symmetric (Z = 0) three terminal junction without bar- 



riers at the NS-interfaces. The continuum current con- 
sists of particles injected from both the normal reservoir 
and the superconducting reservoirs. The current density 
in lead 2 from all injected quasiparticles from the super- 
conductors is 

e 2 sin 8 sinh 7- , , r . 
i s 2 = {sin^ [{AD- 2De- e 2 )cosh7 c 

-2Dee^] - iae^RD - e 2 sin 2 (0/2) cosh 7c } , (Al) 

where Z c = [cos /3(cosh 2-f c (l — e) + esinh27 c ) — R — 
Dcoscj)] 2 + [sin/3(sinh27 c (l — e) + ecosh27 c )] 2 and 7c = 
arccosh(i?/A). In lead 3 we get i^{4>) — — «|( — 4>)- This 
current density is an oscillating function of energy with 
largest amplitude for energies close to E = A and is given 
at negative energies by ij(E) = —ij(—E). 

For the particles injected from the normal reservoir, 
the sum current in lead 2 becomes 



it = 2- 



e e sin (3 cosh 7, 
h Z r 



{sin <p [2D cosh 7c + e sinh j c ] 

+4cr v /i?L> - e 2 /4sin 2 (0/2)sinh 7c | (A2) 

with i^((j>) = — *J(— 4>) in lead 3. The difference current 
in lead 2 has the form 



e e cosh 7 C r 

h Z~ r V 



cosh 7c [e cos (f> + 2a \J RD — e 2 /4 sin < 



+ (1 — e) [sinh 7c — sinh3 7c ] — ecosh3 7c 
-I- cos /3[2 sinh 7c (i? + D cos 4>) 

+ cosh 7c (e(l + cos0) + 2a^jRD - e 2 /4sin(/>)] j (A3) 

with ig(<f>) = — i^(— 0) in lead 3. For negative energies 
we get i+(E) = -ij(-E) and iJ(E) = ij(-E). 



APPENDIX B: SPECTRAL DENSITIES OF 
CURRENTS IN THE WEAK COUPLING LIMIT 

In this appendix we analyze the central quantity in the 
current density expressions (|l6|)-(|l8|), given by 

e e 

Z ~ [(1 - e) cos 26 - Rcos X - D cos <j)} 2 + e 2 sin 2 29 ' 

(Bl) 

in the limit of zero coupling e —> 0. We can conveniently 
rewrite 



1 



(B2) 



with F(E, (j)) = [(l-e)cos26>-.Rcosx-£>cos0]/sm26>. 
In the limit of zero coupling the expression becomes 



lim 



F 2 



nS(F) = J^TT 



d_ 

dE 



F 



8{E-Ei 



(B3) 
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with the energies E^ , given by 

cos 29 — R cos x — D cos = 0, 



(B4) 



being the energies of the bound Andreev states. By 
rewriting 



D sin < 



dE dEdcb dE sin26> 



(B5) 



the expression (|B3|) becomes 

lim — = ■ 
Z A' 



D| sin sin 20 1 



dE 



5{E-Et). (B6) 



From Eq. (B4), the derivative of energy with respect to 
phase becomes 



dE_ 
~d4 



D sin < 



2 sin 26 



( j + + -!— R sin x ^ 

\^/A7^W ~ hv F ~ hv F sin 29 J 



(B7) 



Using the relation (B4) and the fact that R + D = 1 we 
can rewrite 



| sin 26»| = yJ{D + R) 2 - (D cos + R cos X ) 2 



D 2 sin 2 + R 2 sin 2 x + 2DR{1 - cos cos X ) ■ (B8) 



which shows that |sin2#| > i?| sin xl ■ This gives that, 
since X > I by definition, the factor in the parenthesis 
in the denominator in Eq. (B7) is always positive. The 
relation 



'dE . 
sgn [ -r-r sin <p sm 20 



-1 



(B9) 



then follows from Eq. (|B7|). The Eqs. flBJ) and (gg) are 
the technical result of this appendix. 



APPENDIX C: CURRENTS FOR DIFFERENT 
LENGTHS 

Here we list all expressions for the partial currents in 
different length limits. The junction considered is a sym- 
metric (I = 0) three terminal junction without barriers 
at the NS-interfaces, in the weak coupling limit e«l. 



h l-X>sin 2 (0/2) V ' V ' 



ee^jsin(0/2)i 

?il-Xsm 2 (0/2) V ° ; V ; 

When there are two Andreev levels, with 1 -C X>, 
(3/2 > \[T) and X7q" ~ Xq" ~ A, the currents become: 



sA X 



X>sin(0) 



k to 2|sin(0/2)|^/l -£> sin 2 (0/2) 
[tanh^o /2fcT) - tanh(X+ /2feT)] , (C5) 



X = 



;A X 



X> sin(0) 



71 £° 4|sin(0/2)U/l -X) sin 2 (0/2) 



[.g(X+)- 5 (X -)], 
(C6) 



A z!4o 1 — D sm (0/2) 



= C ° Am [H E o) + H E o)}- (C8) 



^ ^° Jl-Xsin 2 (0/2) 



6 _ e ?Wg 

eq ~ h L 



2. Long limit (X < £ ) 

A/Dsin(0) 



AT 



2| sin(0/2)| y/l - X>sin 2 (0/2) 
x ( [tanh(X r //2fcX) - tanh(X+/2fcT)] J 

\n=0 / 

-M*tanh(A/2fcT), (C9) 



( N 

1 - D sin 2 (0/2) \ 71=0 



ehvF vX>cos(0/2) 



h 2L 



g(E-)-g(E+)] 



(CIO) 



1. Short limit (X = 0) r e Twf oVX>i?sm0 v^rwp+x . 

A I " = n — 2 [ l-D S in\ m] Z} h{E « ) + h ^ 

4 b , = V / bm tanh(X -/2fcX), (CI) fc , 

fi 2^/l-Xsin 2 (0/2) ' ' 



e ehvF cos 2 (0/2) v~> r , ^ , ^ 
eA X^sin0 X = - F ^ j X>(£~) + • 

X- = — ■ = =g(E ), (C2) ft 2X 1 - X>sin (0/2) ^ 



4 A /l -X sin 2 (0/2) 

(C12) 
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where h(E) = tanh[(E - eV)/2kT] - tanh[(J5 
eV)/2kT] and g{E) = tanh[(£ + eV)/2kT] + taring 
eV)/2kT] -2tanh(E/2kT). 
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